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ABSTRACT 



We present deep XMM-Newton observations and ESO WFI optical imaging of two X-ray-selected fossil group candidates, 
RXCJ0216. 7-4749 and RXCJ23 15.7-0222. Using the X-ray data, we derive total mass profiles under the hydrostatic equilibrium 
assumption. The central regions of RXC J0216.7-4749 are found to be dominated by an X-ray bright AGN, and although we derive a 
mass profile, uncertainties are large and the constraints are significantly weakened due to the presence of the central source. The total 
mass profile of RXC J23 15.7-0222 is of high quality, being measured in fifteen bins from [0.075-0.75] A500 and containing three data 
points interior to 30 kpc, allowing comprehensive investigation of its properties. We probe several mass models based on the standard 
NEW profile or on the Sersic-like model recently suggested by high-resolution N-body simulations. We find that the addition of a 
stellar component due to the presence of the central galaxy is necessary for a good analytical model fit. In all mass profile models 
fitted, the mass concentration is not especially high compared to non-fossil systems. 

In addition, the modification of the dark matter halo by adiabatic contraction slightly improves the fit. However, our result depends 
critically on the choice of IMF used to convert galaxy luminosity to mass, which leads to a degeneracy between the central slope of 
the dark matter profile and the normalisation of the stellar component. While we argue on the basis of the range of M,/L]t ratios that 
lower M./Lr ratios are preferred on physical grounds and that adiabatic contraction has thus operated in this system, better theoretical 
and observational convergence on this problem is needed to make further progess. 

Key words, galaxy: galaxy group - fossil group: fossil group - X-rays: galaxies: group 



1. Introduction 

. The now well-established cold dark matter (CDM) paradigm lies 

■ at the heart of the fiducial scenario for the formation of structure 
in the Universe. Numerical simulations of hierarchical cluster- 
ing in the cuiTently-favoured ACDM cosmology make a number 
of observationally testable predictions. One example is the ex- 

, istence of a cusped, quasi-universal dark matter density profile 
that is characterised by the scale radius, r^, at which the loga- 
rithmic slope of the pro file is -2, and a dime nsionless concentra- 
tion parameter, c (e.g.. iNavarro et al]ll997h . The latter parame- 
ter exhibits a distinct mass dependence related to the early for- 

' mation of lower- mass haloes in the hierarchical context (e.g., 

' lDuffvetal.ll2()08h . 

Lying at the nodes of cosmic filaments, clusters and groups 
of galaxies are dark matter dominated objects whose prop- 
erties are a sensitive test of these predictions. Indirect evi- 
dence for universality in the underlying dark matter distribution 
in galaxy groups and clusters was indicated by the similarity 
in X-ray surface brightness and temperat ure profiles observed 
with t he ROSAT and ASCA satelUtes (e.g.. iNeumann & ArnaudI 
119991; IVikhhnin et al.lll999l: iMarkevitch et al.lll998l) . However, 
more detailed, quantitative, information can be gleaned from 
examination of the total mass profiles of these systems 



obtained through the hydrostatic equilibrium equation. The 
vastly improved spatial resolution and throughput afforded by 
XMM-Newton and Chandra have now allowed detailed inves- 
tigation of the mass (and thus total density) profiles of mod- 
erately large samples of clusters and groups of galaxie s (e.g., 
Pratt & Arnaud 2002; Pointecouteau et al. 2005; Vikhli nin et al 
2006; Humphrey et al. 2006 
2007 



Buote et al. 2007; Gastaldello et al 
Schmidt & Allen 2007; Sun et al. 2009. and references 
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therein). These and other works have confirmed the cusped, 
quasi-universal nature of the dark matter profile, with a varia- 
tion in concentration roughly in line with predictions. 

However, the presence of a large concentration of baryons 
in the centre of the system can deepen the potential well and 
modify the distribution of dark matter. There has been extensive 
analytical and numerical work on the modification of the mass 
profile of a dark halo induced by the assembly of baryons in 
the inner regions. Early results suggested that simple analytical 
prescriptions based on the conservation of adiabatic invariants 
gave an accurate description of the halo respons e. Follow i ng the 
early work by Barnes & White ( 1984), Blumen thTl et al.l (Il986h 
devised a simple formula to link the dark mass profiles before 
and after the assembly of a galaxy. Given the initial, spherically- 
symmetric enclosed mass profiles of the dark matter, M'^^(t), 

and baryons, /Wj^(r), one may derive the final dark mass profile, 

f f 

M^^{t), once the final baryonic mass profile, Mj^ (r), is specified. 
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The model assumes that dark matter particles move on circular 
orbits before and after the contraction, and that their initial, r, , 
and final, rj , radii are related by the condition; 



(1) 



rfWiirf) + Mdn,] = r,[Md™ + M',(n)], 

where M^/^ = Md^iff) = A#dm('"i) is the dark mass enclosed 
by each dark matter particle. This modification is generally 
termed 'adiabatic contracti on' due to its pa r ameterisation via 
adiabatic invariants (e. g., Eggen et al. 1962', Blu menthal et alj 



1986t iGnedin et all 12004: Duffy etal. 2010). Further work by 
Gnedin et al.l (|2004|) showed, based on hydrodynamic numerical 



simulations, that a modification to using the orbit-averaged ra- 
dius described their results better than the standard circular orbit 
assumption. 

There are few observational tests of this prediction because 
of the difficulty of obtaining accurate mass profiles this deep 
into the core of a given object. Stellar dynamical studies suffer 
from the velocity dispersion anisotropy problem; weak lensing 
cannot measure the mass on these scales (less than 100 kpc), 
and strong lensing arcs are pr eferentially produced in dynam- 
ically disturbed systems (e.g.. iBartelmann & Steinmetjll996h 
which are unsuitable for the detection of the eff'ect. Although 
some progress has re cently been made on all of these fronts, 
(e.g.. lSand et al.ll2008h . in this context. X-ray observations offer 
several advantages conducive to investigation of the mass distri- 
bution in the central regions. In particular, the overall signal to 
noise of a good quality X-ray observation can be significant at 
many hundreds of sigma, and the density squared dependence of 
the emission results in a centrally-concentrated signal ideal for 
investigation of the core regions. The only assumption needed 
is that the intracluster medium (ICM) be in hydrostatic equilib- 
riuni; 

IZappacosta et al.l (l2006h were the first to test for adia- 
batic contraction of the dark matter using X-ray observations. 
Applying several different analytical models (derived from nu- 
merical simulations) to the dark matter profile of the ~ 3.5 keV 
relaxed galaxy clu ster A2589, they f ound no evidence for adi- 
abatic contraction. iGastaldello et alj (|2007|) analysing the mass 
profiles of 16 galaxy groups, found that addition of adiabatic 
contraction did not improve the fit. Building on this result, 
[Humphrev et al. (2006) examined the mass profiles of seven 
early type elliptical galaxies, again finding no strong evidence 
for adiabatic contraction. 

Observational evidence points to the existence of a class 
of dark matter haloes that have the gravitating mass and hot 
gas content of groups, but which are dominated in the op tical 
by a sin gle large early-type g alaxy dPonman et al.l 119941) . As 
noted bv lVikhlinin et alj (Il999h . such 'fossil' systems may have 
been undisturbed for a very long time if they are the end re- 
sult of merging within a normal compact group, and thus they 
may represent the ultimate examples of hydrostatic equilibriurrQ. 
Furthermore, the exceptionally large concentration of baryons in 
the centre of the potential well makes them ideal systems with 
which to search for the effect of adiabatic contraction of the dark 
matter component. 

In the pres ent paper we build on our previous work on 
fossil systems dZibetti et al.l l2009bh by investigating the mass 
profiles of two fossil group candidates, RXC J0216.7-4749 and 
RXCJ23 15.7-0222, using X-ray and optical data. The central 



' Note that some numerical simulations point to the fossil state being 
only a transie nt phase that is terminated by renew ed infall from the local 
environment jvon Benda-Beckmann et an2008h . 



regions of RXC J0216.7-4749 are dominated by emission from 
a bright AGN, and although we recover the total mass profile, 
the constraints are significantly weakened due to the presence 
of the AGN. For RXC J23 15.7-0222, we are able to recover the 
mass profile to high precision. We find that the central regions 
are dominated by the stellar mass component, and using the 
mass profile of the central galaxy from optical data we inves- 
tigate several parameterisations of the total mass distribution. 
We find that the best fitting mass parameterisation is obtained 
when we add a stellar co mponent both in the case of a NFW 
and a iNavarro et al.l (|2004|) profile, and that the addition of adia- 
batic contraction slightly improves the fit. However, the result is 
critically dependent on the assumed IMF used to convert stellar 
luminosity to mass. We discuss the impact of the assumed IMF 
on our results, and argue, based on the range of M^ILr ratios, 
that adiabatic contraction has indeed operated in this system. 

Unless otherwise noted, we adopt a ACDM cosmology with 
Hq = 70 km s-i Mpc-i (i.e., /170 = 1), Qm = 0.3 and Qa = 0.7. 
All uncertainties are quoted at the 68 per cent confidence level. 
In the following, Rs is the radius corresponding to a density con- 
trast of 5 times the critical density at the redshift of the system. 

2. Optical data 

2.1. Observations 

Optical photometry of the t wo systems was c arried out using 
the Wide Field Imager (WFI lBaade et al.|[T999h mounted at the 
Cassegrain focus of the MPG/ESO 2.2 m telescope in La Silla, 
Chile. The WFI is a focal reducer-type mosaic camera which 
consists of 4 X 2 CCD chips, each with 2048 x 4096 pixels and 
a field of view (FoV) of 8.12' x 16.25' (0.238"/pixel). Chips 
are separated by gaps of 23.8" and 14.3" in the R.A. and Dec 
directions, respectively, so the WFI FoV per exposure is equal 
to 34' X 33', with a filling factor of 95.9%. Imaging in the B 
(Ao = 4562.52 A) and R^ (^0 = 6517.25 A, hereafter simply 
R) broad-band filters was performed on September 19th, 2008 
under photometric conditions but with poor seeing (~1.7 and 
1.6" FWHM in B- and R-band, respectively) owing to a strong 
wind. Twelve dithered observations gave a total exposure time 
equal to 1080 s (B) and 720 s (R) and a resultant imaged re- 
gion of 30' X 30'. The two fields of Landolt (1992) photomet- 
ric standard stars SA92 and SA 113 were imaged in the B, V 
(/!() = 5395.62 A), and R broad bands in order to determine 
the photometric zero-point of the night. Twilight-sky flats were 
taken in the evening and following morning for each night of 
observations. 

2.2. Data reduction 

The WFI data were reduced using the data reduction syst em de- 
veloped for the ESO Imaging Survey (EIS Renzini & da Costal 
1 19971) and its associated EIS/MVM image processing library 
version 1.0.1 (Alambic). The data redu ction steps are de- 
scribed in detail in p revious pubhcations (JArnouts et al.ll200ll : 
iMignano et al.ll2007l) . and are omitted here. 

For images of standar ds, source detection an d photometry 
were based on SExtractor (JBertin & Arnoutsll99 6). Magnitudes 
were c alibrated to the Johnson-Cousins system using Landoltl 
(Il992h standard stars whose magnitudes were obtained using 
a 15 arcsec-wide circular aperture, which proved to be ade- 
quate by monitoring the growth curve of all the measured stars. 
Photometric solutions with minimum scatter (the uncertainty is 
equal to 0.03 and 0.05 mag for the zero points in the B and R 
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Fig. 1. i?-band optical images of RXCJ23 15. 7-0222 (left) and RXCJ0216.7-4749 (right), obtained with the Wide Field Imager on the MPG/ESO 
2.2m telescope at La Silla. The red contours start at Icr above background and increase in steps of ^|2, and are derived from the wavelet-smoothed 
X-ray image. 
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Fig. 2. Colour-magnitude diagrams of RXC J0216.7-4749 (top) and 
RXCJ23 15.7-0222 (bottom). In both panels, filled circles denote ob- 
jects within a projected radius of 0.5/?20o; open circles denote objects 
outside this radius. In the bottom panel, red squares denote spectroscop- 
ically confirmed members of RXC J23 15.7-0222, as listed in NED. 



passbands, respectively) were obtained through a two-parameter 
linear fit to 24 photometric points per passband, the extinction 
coefficient being set equal to that listed in the "definitive" solu- 
tion obtained by the 2p2 Telescope Team. In general, zeropoints 
and colour terms are consistent with those obtained by the 2p2 
Telescope Team or by the ESO DEEP Public Survey (DPS) team 
jMignano et al.ll20 07). The final optical images of the two ob- 
jects are shown in Figure [1] 

For the science frames, source detection and photometry 
were based on SExtractor as well. Galaxies can be robustly sep- 
arated from stars down to 7? = 20 from the behaviour of the 
flux radius (i.e., the radius containing 50% of the light) and 
the stellarity index as a function of magnitude. Visual inspec- 



tion is used to test the selection criteria. At fainter magnitudes, 
a conservative separation can be obtained only by assuming a 
threshold value of 0.98 for the stellarity index, below which 
all selected, non-saturated/non-truncated objects are classified as 
bona-fide galaxies. For the brightest and largest galaxies within 
the value of R200 of either gr oup (likely member galaxies of 
RXCJ 23 15.7-0222) the Kron ( lKronll980h magnitudes obtained 
through SExtractor were replaced by the total magnitudes ob- 
tained by analytically extrapolating to infinity the azimuthally 
averaged surface brightness (SB) profiles. These are measured 
by means of elliptical isophote fitting thro ugh the IRAF task 
ellipse as detailed in lGavazzi et al.l (GOOl V Spurious overlap- 
ping stars/galaxies were conservatively masked out before ellip- 
tical isophotes are fitted to the galaxy of interest. Each SB profile 
was eit her fitted using a sing le 'de Vaucouleurs' r''''*-law com- 
ponent (Ide VaucouleursI 19481) or, alternatively, decomposed into 
an exponential law plus an inner 'de Vaucouleurs' r'^'*-law. Fits 
were performed assigning the same weight to each data point 
(logarithmically spaced) in order to properly weight the outer 
parts of the profiles, which would be highly under-weighted in 
a scheme based on the signal-to-noise ratio (S/N). Furthermore, 
points at semi-major axis smaller than the seeing FWHM were 
excluded. Apparent surface brightness profiles and, thus, total 
magnitudes w ere corrected for atm ospheric extinction and galac- 
tic extinction ( ISchlegel et al.|[l99 8). Figure |2] shows the colour- 
magnitude diagrams of the two systems, which are discussed fur- 
ther in Sect. |4] 

For NGC7556, the brightest cluster galaxy (BCG) of 
RXC J23 15.7-0222, the corrected R-band surface brightness 
profile was integrated along the radial distance to give a pro- 
jected R-band luminosity profile, adopting a luminosity distance 
of 104 /i^Q Mpc. This luminosity pr ofile was then depr ojected 
through the algorithm developed by iMagorrianI d 19991) for an 
axysimmetric galaxy, with the kind help of J. Thomas. This al- 
gorithm finds the full range of smooth axisymmetric (luminos- 
ity) density distributions consistent with a given surface bright- 
ness distribution and inclination angle. The best solution corre- 
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spends to the (luminosity) density that max imizes a penalized 
log-likelihood function ('see lMagorrianlfl999l) . 



3. X-ray data 

3.1. Preliminaries 

RXCJ0216.7-4749 and RXCJ23 15.7-0222 were observed by 
XMM-Newton in THIN filter mode for 46 ks and 43 ks, re- 
spectively. Calibrated event files were produced using the stan- 
dard chains in SASvS.O, and data were cleaned for soft proton 
fl ares, PATTERN se lected and vignetting corrected as detailed 
in lPratt et alj (l2007h . The resulting exposure times are 40 ks for 
RXC J0216.7-4749 and 33 ks for RXC J2315.7-0222. 

Point sources were detected on a broad band image via the 
SAS task ewavdetect, with a detection threshold of 5cr. After 
visual checking of the results, these sources were excluded from 
further analysis (excepting the central AGN in RXCJ0216.7- 
4749 which is discussed in more detail in the Appendices). 

3.1.1. Background subtraction 

Once flares are removed, the remaining background can be sep- 
arated into two main components: 

- the particle background, which is dominant at high energy, 
and 

- the cosmic X-ray background (CXB), consisting of, at high 
energy, the unresolved AGN, and at low energy, emission 
from the Local Bubble and the halo of the Galaxy. 

Using a custom stacked background assembled from obser- 
vations taken with the filter wheel closed (FWC), recast to the 
same observation aspect as the source files, we adopt the follow- 
ing procedure for the background subtraction: 

1. We use the count rates in a high energy band ([10-12] keV for 
EMOS, [12-14] keV for EPN) to normalise the FWC prod- 
ucts to those of our observation. 

2. For the surface brightness profiles, we subtract the nor- 
malised FWC data and then determine a region external to 
the group emission where the surface brightness profile is 
flat. From this area we determine the count rate due to the 
CXB, which we then subtract from the surface brightness 
profile. 

3. For the spectral analysis, we extract a spectrum from the re- 
gion external to the group emission, which we fit using a 
double unabsorbed thermal emission model for the galaxy 
and the local bubble, plus a powerlaw with fixed slope 1.4, 
absorbed with the Galactic column density in the direction 
of the group. This background spectrum, scaled by the ra- 
tio of the areas, is then added as an extra component in each 
annular spectrum of the temperature profilq^- 

RXC J23 15.7-0222 fills the XMM-Newton field of view, 
which complicates the background subtraction. We discuss the 
background subtraction for this object in AppendixlAl 

3.2. Gas density distribution 
3.2.1. Morphology 



Figure O shows an XMM-Newton/ESO WFI X-ray/optical 
overlay image for each system. Each fossil group exhibits mor- 



^ Energy bands corresponding to fluorescence lines are ignored dur- 
ing the fit. 



phologically regular emission with the X-ray peak centred on 
the brightest cluster galaxy, as expected if the systems have been 
dynamically quiescent for a considerable period of time. 



3.2.2. Gas density profile 

The gas density profiles were derived as described in 
ICroston et alj (l2008h . Briefly, we extract the surface brightness 
profiles in 3'.'3 bins centred on the peak of the X-ray emission 
using photons from the [0.3-2] keV energy band. After back- 
ground subtraction the profiles are rebinned to a significance of 
3cr per bin. Deprojected, PSF-corrected emission measure pro- 
files were produced u sing two methods: (i) the non-parametric 
method described in ' Croston et al.l (|2006[) and (i i) bv fitting 
a projected 3D parametric model dVikhlinin et al.l 12006). con- 
volved with the XMM-Newton PSF, to the surface brightness 
pro file. In both cases we use redistribution matrices derived from 
the iGhizzardil (l2002h parameterisation of the XMM-Newton PSF. 
Profiles were converted to gas density using a conversion fac- 
tor derived in XSPEC from the global temperature measured 
in the [0.15-l]/?50() aperture. Finally, the temperature and abun- 
dance profiles were fitted with analytic al models (adapted from 
lAllenetaLlbOOU IVikhlinin et"aLll2006h and used to calculate a 
correction factor for each density profile bin due to radial varia- 
tions of these quantities. The surface brightness and gas density 
profiles of each group, derived using the two methods, are shown 
in the left and centre panels of Figures |3] and |4] The parametric 
model is obviously smoother than the non-parametric deprojec- 
tion, but the differences between the results from the two decon- 
volution methods are minimal. 

For RXCJ0216.7-4749, the bright central AGN is evident 
as a point source in the X-ray surface brightness profile, which 
complicates deprojection and introduces significant uncertainty 
into the properties in the very central regions. For this object, we 
first fitted the surface brightness profile with an analytical model 
(Vikhlinin et al. 2006) plus a point source, with the point source 
normalisation constrained from a spectral fit of the central re- 
gion as detailed in Appendix IBJ We then subtract this contribu- 
tion from the total surface brightness profile, adding 5 per cent 
systematic errors to the remainder, and deconvolve as detailed 
above. 



3.3. Spectral analysis 

3.3.1. Spectral fitting 

We extract spectra in circular annuli centered on the peak of X- 
ray emission. Annuli are defined so as to have a signal to noise 
ratio better than 30cr in the [0.3-2] keV energy range. After back- 
ground subtraction they were binned to a minimum of 25 counts 
per bin. We then fitted an absorbed thermal bremsstrahlung 
model plus the background model described in Section [3. 1.1 1 to 
each spectrum using x^ statistics. The background model has 
all parameters fixed and the normalisation of each component 
is scaled to the ratio of the extraction areas. Absorption is fixed 
to the Galactic value in the direction of the group. In the exter- 
nal regions of RXC J0216. 7-4749 the spectra are of insufficient 
quality to reliably constrain the abundances. In these annuli we 
freeze the abundance at 0.3 Z© and fit only the temperature and 
normalisation of the spectra. We fit in the [0.3-10] keV band 
excluding instrumental emission lines. The AGN component in 
RXC J0216.7-4749 is modelled as described in AppendixO 
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Fig. 3. Radial profiles of RXCJ23 15.7-0222. Left panel: Background subtracted, vignetting corrected surfa ce brightness profile 
with best fitting projected analytical model (red) and projected deconvolved profile using the method of dCroston et al.l |20()6|) 
(green) overlaid. Vertical dashed lines illustrate the extent of the annular bins from which spectra were accumulated. Centre panel: 
Corre sponding density profiles. The solid blue line is the density profile derived from the best fitting AB model (Pratt & Arnau^ 
I2OO2I) to the surface brightness profile. Right panel: Temperature profile. The black points with error bars depict the projected tem- 
perature profile. The solid green line is the deprojected, PSF corrected 3D profile; dotted lines show the associated uncertainties. 
The solid red line is the reprojected 3D model. 
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Table 1. Basic data. 



Group RA Dec 


z Lx.bo^ kT 
(ergs"') (keV) 


(xlO'^Mo) 


(kpc) 


{Rsoo) 


RXCJ0216.7-4749 02n6"'45'';0 -47°49'09'.'8 
RXCJ23 15.7-0222 2345"'44^:4 -02°22'54'.'4 


0.064 3.2 X 10"' 2.05 ± 0.05 
0.026 2.1x10*3 1.68 ±0.03 


8.14±0.16 
5.36 ±0.37 


645.7 ± 5.0 
568.7 ± 13.1 


1.0 
0.73 




liming Q.M = 1, i^A = and ho = 0.5. 








" Bolometric ([0.001-100.] keV) luminosity calculated ass 




* Estimated from the M=;rai - Yy relation of|Amaud et alj |2007|). 





3.3.2. Temperature profiles 

The projected temperature profiles of the two groups are shown 
in the right hand panels of Figures [3] and |4l When plotted with 
a logarithmic radial axis, the temperature profiles for both of 
the groups exhibit the bell shape typi cal of cool core clusters 
(IVikhlinin et al.ll2006t IPratt et alJ 12007 ). Deprojection and de- 
convolution of the temperature profile s is undertaken by fitt ing 
parametric 3D models (adapted from IVikhlinin et al] l2006h to 
the projected profiles. These models are convolved with a re- 



sponse matrix that simultaneously takes into account projection 
and PSF redistribution, then projected and fitted to the observed 
temperature profile. I n projecting the models, the weighting 
scheme introduced by ' Vikhliniril (I2006L see also iMazzotta et al.l 
2004) is used to correct for the bias introduced by fitting isother- 
mal models to multi-temperature plasma. Uncertainties are com- 
puted using a Monte Carlo procedure; these are subsequently 
corrected to take into account the fact that parametric models 
tend to over-constrain the 3D profile. Full details of the method 
will appear in a forthcoming paper. 
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4. Reference mass, radius and fossil status 

We make a first estimate of M500 of each system by iteration 
about the M^oi) - Yx relation of'A rnaud et alJ (l2007h . as described 
in|Kravtsov et al. (2006). 7^ is the product of the gas mass inside 
/?5()0 and the temperature in the [0.1 5-0.75] /^soo region, a quan- 
tity which has been shown to be a robust, low scatt er mass proxy 
in numerical siniulatio ns and observations (e.g., lArnaud et alj 
l2007tlPoole et al.l2 007^. The resulting temperatures, masses and 
radii are listed in Table[T] 

The enipirical definition of a fossil group as established by 
aP (l2003h requires a bolometric luminosity L > 10'*^ /i,,? 



E 



ones et 



erg s" in X-rays, and a magnitude gap of Ami2 > 2 mag, where 
Ami2 is the absolute total magnitude gap in R band between the 
brightest and second brightest galaxies in the system within half 
the projected virial radiufl 

The spectroscopic bolometric luminosity of both systems is 
well above the 10^' h^^ erg s"' limit defined by (see Table [TJ. 
Regarding the magnitude gap, the colour-magnitude diagrams 
shown in Figure |2] clearly reveal a large magnitude gap be- 
tween the brightest and second brightest galaxies in both sys- 
tems. Assuming the values for T^soo listed in Table [1] and a ra- 
tio Rioo/Rsm - 1-39 ( Pointecouteau et al. 2005 ). then Ami2 - 
2.21 ± 0.03 mag for RXCJ02 16.7-4749 but is slightly less (at 
< 2cr) for RXCJ23 15.7-0222 (Amn = 1.88 + 0.03D 

However, a system's fossil natur e clearly de pends critically 
on the definition of the virial radius. I Jones et al. used the global 
temperature in combination with a scaling relation derived from 
the numerical simulations in order to calculate their /?2()o; some 
of their global temperatures were estimated from the luminosity 
by assuming a different scaling relation. The R500 (and thus R200) 
measurements we use in this Section are derived from a scaling 
relation derived from X-ray hydrostatic mass measurements, and 
as shown below in Sect. 15.31 M500 varies by 20 per cent or so de- 
pending on the exact modelling of the mass profile, leading to 
changes on the order of an arcminute in the size of the projected 
R200 aperture. Given the above consideratio ns and the fact that 
the fossil criteria fixed bv lJones et al.l (l2003h are somewhat arbi- 
trary, while RXC J23 15.7-0222 does not strictly fulfil the criteria 
within the Icr uncertainties on the photometry, we will continue 
to refer to it as a fossil system in the following. 



5. IVIass profile modelling 

5.1. Mass profile calculation 

We derive the total mass profile of each group from the density 
and temperature profiles (Figures |3] and |4| using the hydrostatic 
equilibrium equation: 



M(<R) 



kT{r)r 



dYaricir) dlnT(r) 

+ 



dlnr 



dlnr 



(2) 



The regularised gas density profiles exhibit structure that is am- 
plified in the radial derivative. To overcome this, we fitted the 
3D gas density profiles with the parametric model described by 
IVikhlinin et alj d2006.) . At each point corresponding to the ef- 
fective radius of the deconvolved temperature profile, we then 
use the radial derivative given by the parametric function fit with 



^ iJones et aTi assume that the virial radius is equivalent to i?200- 
* Spectroscopic redshifts are available in NED for several group 
members in this case and, in particular, for the members brighter than 
R = 14 which fall inside the region imaged with WFI. 



uncertainties on dlogn^/dlogr given by differentiation of the 
regularised density profile. 

Uncertainties on each mass point are calculat ed using a 
Monte Carlo approach based on that of iPratt & Arna ud (2003|). 
A random temperature is generated at each radius at which the 
temperature profile is measured, and a cubic spline is used to 
compute the derivative. In the randomisation, we only keep pro- 
files that are physical: they must increase monotonically with 
radius and the randomised temperature profiles must be convec- 
tively stable, i.e., dlnT/dlnn,, < 2/3. In total 1000 such Monte 
Carlo realisations were calculated; the error on the derivative is 
then the region containing 68 per cent of the realisations on each 
side. The resulting mass profiles are shown in Figure|5] 

5.2. Mass profile models 

We investigated various parameterisations of the total mass pro- 
file, as listed below. 

NFW : We first fitte d these profiles with the integrated mass 
from an NFW profile dNavarro et al.lll997h . viz: 



M(< r) = 4npciz) 6,- r^ 



ln(l+r/r,)- 



r/r, 
1 + r/rs 



(3) 



where pdz) is the critical density of the universe at redshift z, 
Ks is the scale radius where the logarithmic slope of the den- 
sity profile reaches -2, and 6c is a characteristic dimensionless 
density. This model has been shown to be an adequate fit to the 
mass profiles of many morphologically relaxed systems above 
~ 2 keV (e.g.. 'Pra tt & ArnaudlllOOl iPointecouteau et al.ll2005l: 
Vikhlinin et ahaOQ^. 



N04 : We also considered the universal profi l e deri ved from 
higher resolution simulations by iNavarro et alJ ('2004'). The in- 
tegrated mass of this profile is given by (e.g., IZappacosta et al.l 
120061) : 

M{< r) = 4Kp^2r'_2 ^ ^^pf^) (f f "^(^' ^"^"j ^4) 



where 7(77, /I) - J. t'' ^e ' dt is the lower incomplete gamma 
function, and x - r/rs. As discussed below in Section 15.3.21 we 
fit both w ith a free and hmited to lie within the + Icr uncertainties 
found bv lNTvarro et all (l2004 (i.e., 0.14 < a < 0.20). 

NFW+star and N04+star : Previous investigations of group- 
scale haloes have suggested that below a certain mass scale (typ- 
ically ~ 2 keV), the BCG begins to make a substantial contri- 
bution to the distributi o n of mass in th e central regions (e.g., 
iGastaldello et alJ l2007t ISun et alJ l2009l) . As discussed in the 
Introduction, this contribution is potentially even more impor- 
tant in the case of a fossil group, where there is a considerable 
concentration of baryons in the centre of the potential well. We 
thus fitted the total mass profile with the NFW and N04 mod- 
els detailed above with the addition of the mass from the cen- 
tral galaxy derived from our optical observations as described 
in Sect. |2] The fits are initially undertaken with a fixed normal- 
isation for the stellar component of M^ILr - 1.84, as obtained 



V 




J. Democles et al.: Adiabatic contraction in fossils 

10° 



10" 



K 10 
V 

10" 



-2 




^ — I I I I 



^ — I — I — I I I 'I 1 1 1 — I — M-t- 



f--^-^-^-f-i-iUi-f-f-i-^-^| 



I 



10 100 

Radius (kpc) 



10 100 

Radius (kpc) 



Fig. 5. Hydrostatic total mass profiles for RXCJ2 3 15.7-0222 (left) and RXC J0216.7-4749 (right). Each profile is renormalised to 
the mass expected from the M500 - Yx relation of lArnaud et al. (2007) . Vertical dashed lines show r = 0.15 R500 and r - Rsoo- The 
solid red Une is the best fitting NFW mass profile in each case. 



from the B - R colour and the formulae given in IZibetti et alJ 
(2009a), based on a Chabrier (2003) IMlQ. 

In these fits, we first subtract the gas mass from the total 
gravitational mass profile to isolate the dark matter and stellar 
components. We then fit the resulting dark matter plus stellar 
mass with an NFW profile plus the stellar mass profile derived 
from our optical data (Section |2]). Since we are fitting the dark 
matter plus stellar mass profile, the concentration we derive is 
for these components only. In order to calculate the total mass 
concentration including the mass of the hot gas, we add the mass 
of th e hot gas to that of the d ark matter plus stars, and iterate (see 
e.g.. lZappacosta et alj|2006l) . 



NFW+A: star and N04+^ star : Adoption of different IMFs may 
result i n mass estimates th at differ by up to a factor 2, as shown 
by e.g.. lBruzual & Charloti (12003.. their Figure 4). The exact form 
of the IMF is still a matter of debate (see Section 16.1.21 ) and 
therefore we also investigated the influence of freeing the nor- 
malisation of the stellar mass component. The total mass and 
corresponding NFW parameters are estimated by iteration as de- 
scribed above. We discuss the impact of the assumed IMF in de- 
tail in Sect. |6T2l 



NFW*AC+/tstar and N04*AC+A:star : We also investigated 
the possibility of adiabatic contraction o f the dark ma tt er pro - 
file using the prescription described in iGnedin et alJ (l2004l) . 



^ By directly reverting to the model library of IZibetti et al.l ( l2009al) 
and excluding models with amounts of dust that are unrealistic for early- 
type galaxies, we obtain a median likelihood M,/Lf. = 1.41 (thus 0.43 
times lower than obtained with the above-mentioned fitting formula) 
with a ±1 cr(3 cr) range of 1.03-1.80 (0.76-2.27). The adopted value 
of 1.84 should thus be regarded as an upper limit for a Chabrier IMF. 
Adopting a Salpeter (1955) IMF, all M/L ratios would be scaled up by 
a factor of 1.75. 



as implemented in the adiabatic contraction code contrcQ 
We c h oose this p r escrip tion among others ( Blumenthal et al] 
Il986l : lAbadi et all l2009h for consis te ncy with previous in- 
vestigati ons (Zappacostaet al] 120061 : iHumphrev et al.l 120061 : 
Gastalde llo et alJl2007h . In addition, contra has been shown to 
perform well, in a statistical sense, in simulations of the dark 
matter response to galaxy formation for galaxies with total stel- 
lar masses at z = in the range 2.5-6 x lO'" h M^ havi ng ex- 
perienced little star formation over the last 8 Gyr dTissera et al.l 
2009). 

For input to the contra code, we need the initial and final 
baryon profile and the initial dark matter profile. In this process 
only the stellar component is considered to affect the dark mat- 
ter profile. We then proceed as follows: from the observed total 
mass profile we subtract the observed gas mass profile. We then 
fit with an NFW+fcstar or N04+fcstar model. Further subtract- 
ing the stellar mass component after this fit gives the initial dark 
matter profile used as input to contra. The initial baryon (stel- 
lar) profile is assumed to follow the same form as the initial dark 
matter profile, with a normalisation given by: 



(5) 



A^stellar,imt(< '") = /stellar X MDM,init(< '") 

where /steiiar is the stellar fraction at the maximum radius of the 
optical data. The final baryon profile is given by the observed 
stellar mass profile with appropriate normalisation, contra is 
then run on the dark matter profile to compute the adiabatic con- 
traction due to the stellar mass component. The resulting total 
mass profile is then compared to the data and the model iterated 
until convergence. 

As we will see below, the major contributions to x^ come 
from the inner data points (R < 30 kpc), where the uncertainties 
on the temperature and density are smallest. Many of the mod- 
els are formally acceptable, and it is not possible to distinguish 
between them purely in terms of ;i^^. FoUowing the discussion 



ihttp : //www ■ astronomy . ohio- state . edu/~ognedin/contra/ | 
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Fig. 6. Additional modelling of the mass profile of RXCJ23 15.7-0222. Left panel: NFW model fitting. The solid line is the best 
fitting NFW model. The dashed line adds the stellar component with free normalisation (NFW+A:star); the dot-dot-dot-dashed 
line adds adiabatic contraction to the dark matter component (NFW*ACH-A:star). Right panel: N04 fits, with same line styles. 
Corresponding best fitting parameters are listed in Table|2] the quality of the fit are roughly the equivalent with and without adiabatic 
contraction; the main differencies are in the normalisation of the stellar mass profile. 

Table 2. Results of fits to the mass profile of RXCJ23 15.7-0222. 



Model 



Msoo (10'-^ Mo) R500 (kpc) 



/gfli.500 



MJLr MioolM, 



ifflL. 



X Idof dvin 



NFW 


8.01 ± 0.82 




3.58 ±0.29 


497 ± 


14 


0.093 ±0.013 




0.67 ± 0.07 


17.9/13 


1.08 


NFW+star 


7.33 ±0.82 




3.99 ±0.33 


515± 


14 


0.089 ±0.013 


1.84 


0.74 ± 0.08 


7.9/13 


0.72 


NFW+istar 


5.68 ± 0.93 




4.33 ±0.37 


530 ± 


15 


0.086 ±0.013 


4.40 ±1.08 


0.81 ±0.09 


2.8/12 


0.14 


NFW*AC+star 


6.01 ±0.59 




3.99 ±0.27 


515 ± 


12 


0.089 ± 0.012 


1.84 


0.74 ± 0.07 


3.65/13 


0.15 


NFW*AC+ k star 
a free 


5.61 ±0.12 




4.17 ±0.27 


523 ± 


11 


0.087 ±0.011 


1.81 ±0.33 


0.78 ± 0.07 


2.67/12 


0.14 


N04 


5.87 ±0.97 


0.108 ±0.016 


4.30 ±0.35 


528 ± 


14 


0.085 ± 0.012 




0.80 ± 0.09 


4.63/12 


0.22 


N04+star 


5.53 ± 1.37 


0.138 ±0.042 


4.42 ±0.31 


533 ± 


13 


0.085 ±0.011 


1.84 


0.82 ± 0.08 


3.78/12 


0.15 


N04+j(: star 


5.30 ± 1.04 


0.299 ±0.153 


4.03 ±0.31 


517 ± 


14 


0.088 ±0.012 


4.50 ±2.34 


0.75 ± 0.08 


2.36/11 


0.14 


N04*AC+star 


5.07 ± 0.95 


0.241 ± 0.074 


4.20 ± 0.30 


524 ± 


13 


0.086 ±0.012 


1.84 


0.78 ± 0.08 


2.53/12 


0.14 


N04*AC+i star 


5.12±1.14 


0.326 ±0.182 


4.05 ± 0.30 


518 ± 


13 


0.088 ± 0.012 


2.09 ±1.18 


0.76 ± 0.08 


2.09/11 


0.13 


0.14 < Q- < 0.20 limited 






















N04 


6.81 ±0.82 


0.142 ±0.005 


4.01 ±0.31 


516± 


13 


0.088 ±0.012 




0.75 ± 0.08 


6.78/12 


0.43 


N04+star 


6.12 ±0.85 


0.154 ±0.021 


4.28 ± 0.26 


527 ± 


11 


0.086 ±0.011 


1.84 


0.80 ± 0.08 


4.02/12 


0.16 


N04+/: star 


5.39 ± 1.05 


0.189 ±0.025 


4.30 ±0.28 


528 ± 


12 


0.086 ±0.011 


3.64 ±1.39 


0.80 ± 0.08 


3.05/11 


0.14 


N04*AC+star 


4.83 ± 0.65 


0.193 ±0.021 


4.43 ± 0.25 


533 ± 


10 


0.084 ±0.010 


1.84 


0.83 ± 0.08 


3.29/12 


0.13 


N04*AC+j(: star 


5.32 ± 1.19 


0.191 ±0.024 


4.35 ± 0.26 


530 ± 


11 


0.085 ±0.011 


1.51 ±0.71 


0.81 ± 0.08 


3.01/11 


0.14 



in iGastaldello et al.l (l2007h . we thus also give the values of the 
maximum fractional deviation dvi^ax, which gives equal weight 
to all radial bins, as an additional discriminator of the goodness 
of fit. The dvimax was extensively used by iJin^ (I2OOOI) in an in- 
vestigation of the goodness of fit of the NFW analytical form to 
dark matter haloes both in and out of equilibrium, where it was 
found that a dvimax < 0.3 represented a good fit of the NFW 
model. 



5.3. Mass fitting results 



5.3.1. RXCJ021 6.7-4749 



While the weighted effective radius jLewis et al] 120031 ) of the 
outermost temperature point lies at ~ R500, the uncertainties on 
the mass profile of RXCJ0216. 7-4749 are severely affected by 
the presence of the central AGN. As a result, a simple NFW 
model with C500 = 2.22 ± 0.25 and M500 = l-31!o]i x IO'^^Mq 
is an acceptable fit (^^/dof = 9.8/13) to the mass profile of 
this object once corrected for the AGN contribution. The AGN- 
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corrected profile and best-fitting simple NFW model are shown 
in the right hand panel of Figure |5] 

Partly the good fit is because of the large uncertainties in the 
inner regions which are the result of the correction for the central 
AGN. Appendix ICJ compares the parameters of the NFW model 
fit to the mass profile of RXC J02 16.7-4749 in the cases of (i) no 
correction for the AGN contribution; (ii) excision of the central 
18", and (iii) modelling of the AGN contribution as described in 
Appendix|B] The total mass varies by 30 per cent, with the low- 
est value derived for the uncorrected mass profile and the highest 
value derived for the mass profile corrected for the AGN con- 
tribution. The concentration parameter, unsurprisingly, varies in 
the opposite sense, being highest for the uncorrected profile and 
lowest for the corrected profile. 

The fact that we detect the mass profile directly to ~ R^m 
gives us confidence that the result is not due to an extrapolation 
bias resulting from fitting the mass profile in a limited radial 
range. However, we do not pursue more complicated model fits 
to the mass profile of this system due to the limited quality of 
the data within 0.15 ^soo- 

5.3.2. RXC J231 5.7-0222 

In contrast, the mass profile data of RXC J23 15.7-0222 are of 
sufficiently good quality for a more detailed investigation: there 
are three well-constrained data points within 30 kpc, and the data 
extend to ~ 0.8 /?50(), allowing us to constrain both the mass dis- 
tribution in the centre and the slope of the mass profile in the 
outer regions with good accuracy. Figure |6] shows the mass pro- 
file data and the various model fits which we will now discuss. 

The mass profileof RXC J23 15.7-0222 is not well described 
by the simple NFW model: the reduced ;^f^ is rather high (^'^/dof 
= 17.9/13), and there is a systematic divergence of the data 
points from the model in the inner ~ 100 kpc (Figures |5] and |6]). 
The excess of mass in the inner regions with respect to the NFW 
profile is significant (~ 2.5cr for the inner point), suggesting that 
there is a substantial contribution in the central regions due to the 
stellar mass component. Indeed, the initial fit of an NFW+star 
profile with fixed normalisation of the stellar component results 
in a substantial improvement in x^ relative to the simple NFW 
fit, and the central regions in particular are noticeably better 
fitted. The addition of the stellar mass component results in a 
decrease in the value of the concentration parameter, although 
this is not significant within the relatively large uncertainties. 
Further freeing the normalisation of the stellar mass component 
(NFW -t-A: star) results in another improvement in ;i'^, which an 
f-test confirms is very highly significant (probability 5.3 x 10""*). 

The N04 mass profile allows more centrally concentrated 
mass distributions than the original NFW parameterisation. 
Fitting this form to the total mass data with the a parameter 
left free results in a substantial improvement in the fit statistic 
over the simple NFW case (;^^/dof - 4.63/12), due primarily to 
a better fit of the inner three data points. However, the inferred 
value of a - 0.108 + 0.016 is sig nificantly lowe r than t he mean 
value of 0.172 + 0.032 found bv iNavarro etaP (|2004 . imply- 
ing a density profile that is steeper in the centre and shallower 
at large radius than that inferred fro m CDM. Limit i ng a to lie 
within the ± 1 cr uncertainties found bv iNavarro et alJ (12004 ) (i.e., 
0.14 < a < 0.20) results in a marginally degraded fit (;(f^/dof - 
6.78/12). 

An N04+star model with a free and a fixed stellar normal- 
isation results in a further improvement in x^ over the simple 
N04 case and more interestingly, this fit yields a value of a that 
is in better agreement with the predictions from CDM. Further 



freeing the normalisation of the stellar component results in an- 
other improvement in x^ but at the expense of a higher value of 
a and a higher M^ILr, although then neither parameter is well 
constrained. A similar trend is found f or fits with a Umited t o lie 
within the ± Icr uncertainties found bv iNavarro et alj (|2004|) . 

Finally, we investigated the case where adiabatic contrac- 
tion of the dark matter is taken into account. In both the NFW 
and N04 cases, adiabatic contraction of the dark matter profile 
improves the fit slightly relative to that without. Formally the 
best fitting model is an NFW profile with free stellar normali- 
sation and adiabatic contraction applied to the dark matter, but 
the equivalent N04 profile is statistically almost as good a fit. 
We discuss the implications of this result, and in particular its 
dependence on the assumed IMF, in the next Section. 



6. Discussion 

6.1. Total mass profile of RXC J23 1 5. 7-0222 
6.1.1. Introduction 



In Section |5. 3. 2l above we showed that the best fitting mass pro- 
file model for RXC J23 15.7-0222 requires a central stellar mass 
component, that the fit improves when the normalisation of this 
component is left free, and that the addition of adiabatic contrac- 
tion improves the fit slightly. We emphasise that the radial reach 
of our observations is sufficient to place strong constraints on 
the overall form of the mass profile from deep in the central re- 
gions (< 10 kpc), out to a significant fraction of R^oq. However, 
the main conclusions are driven by the fit to the data interior to 
~ 30 kpc (i.e., the inner four data points), and the effects under 
discussion (i.e., the normalisation of the stellar mass component 
and the effect of adiabatic contraction on the dark matter pro- 
file) are somewhat degenerate. In the following, we attempt to 
disentangle the influence of each effect. 

Given that additional stellar mass is clearly needed, the ques- 
tion of the normalisation of the stellar component becomes cru- 
cial in interpreting the inner dark matter density slope. While 
initially the normalisation of the stellar mass component was set 
to MJLr - 1.84 as detailed in Sect.|521 formally the best fits 
(both in terms of x^ and in terms of dvi,nax) are given by the 
NFW* AC -I- k star and N04*AC + k stai" models where the stellai" 
normalisation is left free and there is adiabatic contraction of the 
dark matter. 



6.1 .2. On the choice of IMF and its impact on the best fitting 
mass model 

As Table |2] shows, the NFW and N04 mass profile fits to 
RXC J23 15.7-0222 significantly improve when stars are in- 
cluded but no prior is put on the IMF. Fits are better in an 
absolute sen se with adiabatic contraction as implemented in 
iGnedin et aLf s code; however, if there is no adiabatic contrac- 
tion, the 3cr ranges for M*/Lr (and thus the IMF) lie between 
~ 1.2 and ~ 7.7 (for the NFW profile) or and ~ 11.5 (N04 
profile), i.e., it is largely unconstrained. This is to be compared 
with the ±3cr confidence range for M^/Lr derived f rom the ob 
served B-R color using the large suite of models by Zibetti et al.l 



(2009a). The predicted range is of 0.76-2.27 if a'Chabrier( 20031 
IMF is assumed, while it is 1.33-3.97 for a'Salpeter ( 1955) IMF 
Thus, with a central value of M*/Lr x 4, the models without 
adiabatic c ontraction seem to favour more bottom-heavy IMFs 
like that of lSalpeteil(ll955h . 
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Interestingly, when an NFW profile is assumed and adia- 
batic contraction is applied to the dark matter, the uncertainty 
on Mi, /Lji drops by a factor of three, and the robust best-fitting 
value of M-i,/Ljt, is remark ably consistent with the value pre- 
dicted assuming aJ Chabrien (12003,) IMF. Similar conclusions can 
be drawn in the case of an N04 profile, but they are much weaker 
since the uncertainties on the values of M* /Lr become larger In 
other words, inclusion of adiabatic contraction leads naturally to 
an M*/Lr that is consistent with a Chabrier IMF whereas the 
models without adiabatic contraction tend to prefer a Salpeter 
IMF 

Hence, the assumption of a particular dark matter profile 
seems to impact on the robustness of the conclusions on the 
IMF. Distinguishing between IMFs at the low-mass end is an 
extremely challenging task: in fact, although low mass stars 
contribute significantly to the mass, their optical/near-IR flux is 
negligible and thus photometric measurements at these wave- 
lengths provide very weak constraints. While the inclusion of 
more bands (especially near-IR) would partly restrict the possi- 
ble range of M^ILr at fixed IMF, this would not alleviate the 
systematic uncert ainty der i ving f rom the unknown IMF. 

The standard ISalpeterl (|1955|) IMF corresponds to a single- 
slope power law ^{M) oc M~* for 0.1 < M < 100 M©, where s - 
2.35. It provides the lowest value of the stellar mass-to-B-band 
luminosity ratio (M-i,/Lb) with respect to analogous IMFs with a 
very steep (dwarf dominated) or a very flat (remnant dominated) 
slopes. Yet for nearby elliptical galaxies, suitable simple stellar 
population (SSP) models with ages of 12 Gyr predict a value 
of M*/Lb which is twice as large as that inferred from dynami- 
cal models (un der the assumption of a constant stellar mass-to- 
light ratio, e.g.. Ivan derMare]||1991h . This excludes any single- 
slope IMF and enforces a flattening of the IMF with respect t o 
flie Salpeter slope below ~ 0.5-0.7 M© (e.g., iRenzinil 120051) . 
Consistently, direct stellar c ounts in Galactic globular clusters 
(iParesce & De Marchill2000l) an d young clu sters with ages rang- 
ing from a few My r to 1 Gyr dde Marchi et al...2005; see also 
iBouvier et alJl2005h point to an IMF with a log-normal form be- 
low 1 M0. This conclusion holds for the present day mas s func- 
tion of fields in the Galactic Disc (Chabrier 2003; Morau xet alj 
12004 ) or Bulge (IZoccal i 2005, and references therein). 

Interestingly, in an analogous Chandra study of mass pro- 
files in seven ell iptical galaxies with either galaxy-scale or 
group-scale halos dHumphrev et al.ll20 06'). the stellar mass-to-K- 
band luminosity ratio (M*/Lk) was_found to be consistent with 
SSP models assuming a lKroupal(l2001 ) IMF. T he latter contains 
a flattening below ~ 0.5-0.7 M©, similar to a IChabried (l2003h 
IMF 

The IMF is also indicated as the most significant source 
of systematic uncertainty in a recent test of adiabatic contrac- 
tion using profiles o f 75,086 elliptical galaxies from the Sloan 
Digital Sky Survey ("S chulz et al.ll200 9). This study is based on 
weak lensing observations in the outskirts of the halo and mea- 
surements of the stellar velocity dispersion in the inner regions 
of galaxies for stacked galaxy samples. Schulz et al. conclude 
that stellar masses need to be larger by a factor of two with re- 
spect to those obtained with a Rroupa IMF to explain the inner 
dynamical-mass excess in their data without adiabatic contrac- 
tion, but su ch an increase would c reate tension with results from 
SAURON (ICaDDe llarietani2006l) . 

From this discussion, we conclude that fits to the mass profile 
of RXCJ23 15.7-0222 yielding low values of M*/Lr (i.e., with 
a Chabrier-like IMF) must be preferred on physical grounds. 
This implies that some form of adiabatic contraction has to be 
invoked whatever the underlying dark matter profile. This con- 
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elusion is at odds with that of Humph rey et al] (12006 *). who cast 
some doubt on the .Gnedin et aL (.20041) adiabatic contraction 
scenario since their best-fitting NFW*AC-i-star models to early 
type galaxy mass profiles yielded significantly lower values of 
M*/Lk than predicted by a Kroupa (2001) IMF. At variance with 
their best-fits, we do see a significant change in the uncertainties 
associated with the best-fitting values of M* /Lr when adiabatic 
contraction is applied, at no cost of tension with the values pre- 
dicted by a range of IMFs. However, we do confirm that allowing 
adiabatic contraction does not produce evident improvements in 
the significance of the best fits, irrespective of the halo profile. 

The different conclusion regarding the role played by adia- 
batic contraction may result from the limited sp atial extent of the 
halo and galaxy regions probed in the study of Humphrey et al] 
(^2006); alternatively, the dissimilar nat ure of the objects under 
study in the present work and in that of [Humphrey et al] (early- 
type galaxies vs. fossil groups) may cause other conclusions to 
be drawn. 

6.1 .3. Characteristics of tlie dark matter profile 

As mentioned in the introduction, fossil systems are so named 
because they are supposed to be the endpoint of the merger his- 
tory of an early-forming compact group and as such their dark 
matter profile should reflect their early age of formation in a 
higher concentration than average. An important question that 
remains to be addressed is whether the concentration of the dark 
matter profile of RXC J23 15.7-0222 is any different from typical 
literature values for non-fossil systems of a similar mass. 

The simple NFW fit to the total mass density profile leads 
to the highest con centration of any of the m odel fits, supporting 
the conclusions of iMamon & Lokai (l2005h . who suggested that 
neglect of the central stellar mass in the central regions, together 
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with a restricted radial fitting range, could mimic a high mass 
density concentration, a s was found by some earlier work (e.g., 
iKhosroshahi et al.l 12004 fl We find that depending on the mass 
model, the resulting values of M500 only cover a relatively small 
range from 3.6-4.4x 10'^ Mq, with corresponding concentration 
parameters c = 5 - 8. Applying adiabatic contraction to the dark 
matter generally lowers the concentration (see Table|2]i. 

Figure |7] shows the concentration parameter derived from 
the various NFW model fits compared to the data compilation 
of relaxed systems from iBuote et aP (1200 7 rl The two fossil 
groups in the present analysis do not exhibit particularly high or 
low concentrations in view of the considerable dispersion in the 
range of measured values. Furthermore, the small span of values 
we find for both c and M500 suggest that we have suflicient radial 
leverage to make an unbiased estimate of these parameters. 

6.1 .4. Effect of possible central point source 

The mass profile results discussed above for RXC J23 15.7-0222 
were derived assuming that there is no additional central point 
source that emits in X-rays. As this might actually be the case, 
we investigated the possibility that a second component may ex- 
ist by fitting the spectrum of the inner three annuli with a power 
law of fixed slope V - 1 .4 in addition to the thermal emission 
model. In the innermost spectral region, a better fit is obtained 
with the MeKaL plus powerlaw model (confirmed by an F-test), 
suggesting that a central X-ray source may indeed be present 
and contributing < 10 per cent of the emission in that regiorQ. 



^ Note that the ext r emely high concentration found for NGC 6482 
by IKhosroshahi et al] ( 120041) may also have exacerbated by the limited 
radial range of their data. 

>* We compute A,,,,, using A,(z) = 18; r^ + 82x[n„(z)-l] -39x[0„,(z)- 
\f with n,„(z) = n,„(z = 0) = 0.3 fsee lBrvan & Normaijrr998') . 

' The upper limit to the X-ray liminosity of this source in the 2-10 
keV band is L ~ 2.9 x 10^' erg s"' which is consistent with the lumi- 



We thus corrected for the effect of this central source on the den- 
sity and temperature profiles as described in Appendix |B] and 
refitted the mass profile with the mass models described above 
in Section 15.21 The best fitting model parameters are given in 
TablelDJi 

While there is a slight change in concentration towards lower 
values (as expected), the most important result is that the same 
trends and conclusions are valid for this analysis as for the anal- 
ysis assuming no central source. We are thus confident that our 
conclusions regarding the properties of the mass and dark mat- 
ter profiles of this system are robust to the presence of a central 
X-ray point source (if any). 

6.2. Gas properties 

6.2.1 . Scaled temperature and abundance profiles 

The scaled temperature and abundance profiles of the two fossil 
groups are shown in Figure [8] When plotted with a logarithmic 
radial axis, the temperature profiles for both of th e groups exhibit 
the bell s hape typical o f cool core clusters (e.g., IVikhlinin et al.l 
2006; Pr att et al.ll2007h . The temperature peak is found interior 
to 0.1 R50Q in each case, further in th an is observed i n clusters 
and in agreement with the findings of Sun et aP (l2009l) from the 
analysis of a large sample of groups. Outside the core regions 
the temperature profiles both decline with a similar slope. We 
fitted the combined data set in the [0.15 - l]/?5oo radial range 
with the model T/T2500 = A^ x (1 -H r/T^soo)", for which we 
find A^ = 1.33 ± 0.03 and a = - 1.13 ± 0.05. Th ese values 
are consistent with those found by ISun et al] (l2009h . suggest- 
ing that fossils are not exceptional in the properties of their in- 
ner or outer temperature profiles when compared to other sys- 
tems. The abundance profiles of the two groups are substantially 



nosity expected from low mass X-ray binaries according to lGrimm et al.l 
12003) study of X-ray binaries in the Galaxy. 
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similar, exhibiting an off'-centre abundance peak similar to that 
seen in some other observa tions of group scale objects (e.g., 
iRasmussen & Ponmanll2007h . where the mean profile rises to- 
wards the central regions, then exhibits a drop in the inner ~ 20 
kpc. 

6.2.2. Entropy and gas mass fraction 

The entropy profiles of the two groups are shown in the left 
hand panel of Figure |9l As is conventional, 'entropy' is calcu- 
lated from the density and temperature profiles: K - kT/n^'^ ; in 
addition, they have been scaled by the characteristic entropy of 
the system 



K,, 



106 keV cm" 



M,, 



lO'- 



Ihi Me 
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E(z)-"' h 



2/3 , -4/3 



(6) 



(e.g.. lVoit et ai]|2005l: iPratt et alj|2010h . The dashed fine shows 
the best fitting power law fit to the median entropy profile in the 
radial range [0.1 - 1 ] Riga for the clust ers formed in the non- 
radiative simulations of IVoit et al] (l2005h . The observed profiles 
exhibit much the same form, and lie very significantly above the 
baseline prediction, as is expected if non-gravitational processes 
affect the ICM to a greater extent at the group scale. The en- 
tropy excess is significant even at large radius (and in the case of 
RXCJ02 1 6.7-4749, this is true at ^500), in agreement with the 
results of I Sun et al] (l2009h and in contrast t o the relative lack 
of excess at lar ge radius fo und in clusters (iNagai et al.l l2007t 
ISun et alJ I2OO9I: iPratt et all l201d ). The radial entropy slope is 
shallower than the typical value of 1 . 1 except in the very outer 
regions (r > Q.6Rsm)- 

The middle panel of Figure |9] shows the integrated gas mass 
fraction profiles /ga,s(< R) - M^^^< R)/M(< R). The gas mass 
fraction increases with radius in each case. The average gas mass 
fraction at 7?2500 is /ga,s,2500 = 0.045 + 0.005, consi stent with 
the res ults of the larger sample of groups studied by ISun et all 
(l2009h . As Table 12] shows, measurement of the gas mass frac- 
tion of RXCJ23 15.7-0222 at /^soo is somewhat model depen- 
dent. The simple NFW model, which gives the worst fit to the 
overall mass profile, yields the highest value of /gas.soo due to 
the model's systematic underestimate of the data point at ~ 400 
kpc (Figure |6l). The most reliable values can be obtained from 
the best fitting models, yielding fgas.so n ~ 0.087 for th is group, 
a value consistent with those found bv lSun et alj (l2009l) for sys- 
tems of similar temperature. 



The link between increased entropy and the to tal gas content 
was recently demonstrated by iPratt et al] (l2010l) . who showed 
that multiplication of the scaled entropy profile with the scaled 
gas mass fraction profile, in effect correcting the entropy for the 
difference in total gas content with radius, yielded entropy dis- 
tributions that were in good agreement with the predictions from 
adiabatic simulations. We show the corresponding fossil group 
profiles in the right hand panel of Figure |9j clearly, this correc- 
tion also works on the group scale, providing further evidence 
that gas content is the key to understanding the physical pro- 
cesses responsible for the behaviour of the entropy. 



7. Conclusions 

We have presented the first deep X-ray and wide field op- 
tical imaging observations of two candidate fossil groups, 
RXC J0216.7-474 9 and RXCJ23 15.7 -0222. Based on the crite- 
ria established by I Jones et aP (l2003h . and taking into account 
the various uncertainties involved in the definition of /?50(), we 
argue that RXC J0216. 7-4749 is a bona fide fossil system and 
RXC J23 15.7-0222, if not formally fossil, shares strong physi- 
cal similarities with this type of object. The X-ray data quality 
is exceptional for this type of object, extending from [0.01 - 
0.75] RsoQ in both cases, allowing us to investigate in detail the 
properties of their profiles. While unfortunately the central re- 
gions of RXC J0216.7-4749 are contaminated by a bright X-ray 
point source that contributes ~ 40 per cent of the emission in the 
central temperature profile bin, we devise a method to correct 
for its presence and derive the resulting corrected density and 
temperature profiles. 

The object temperatures are 2.05 + 0.05 keV and 1.68 ± 
0.03 keV for RXC J0216.7-4749 and RXC J23 15.7-0222, re- 
spectively, when measured in the [0.15-1] /?5oo region, placing 
them squarely in the galaxy group category. Both systems ex- 
hibit regular, highly peaked X-ray emission centred on the BCG, 
indicative of their being morphologically relaxed objects. The 
temperature profiles both describe the shape typical of cool core 
systems, but with a temperature peak at ~ 0.075 Rsm, closer to 
the centre than is observed for more massive systems. Their en- 
tropy profiles show a considerable excess above the expectations 
from non-radiative simulations across the entire measured radial 
range (this is true out to R50Q for RXC J0216.7-4749). 

Using the temperature and density profiles, and assuming hy- 
drostatic equilibrium, we calculated the total mass profiles of the 
two groups. For RXC J0216.7-4749, the best fitting NFW model 
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yields cjoo = 2.22 ± 0.25 and M500 = l-31+°[^ x 10'^ ^^. how- 
ever, mass constraints overall are weak due to the uncertainties 
associated with correction for the AGN, and we do not fit more 
complex models to these data. 

The mass profile of RXC J23 15.7-0222 is of sufficient qual- 
ity for deeper investigation. We find that consideration of the 
stellar mass of the central galaxy is essential to provide a good 
fit to the data. The best fitting m ass model is either the Sersic- 
like profile proposed bv lNavarro et al. (2004) with an index a in 
agreement with predictions, or an NFW profile, in each case plus 
a stellar component. The concentration is not especially high 
compared to non-fossil systems, and appears to be in the range 
observed and expected for normal systems. Applying adiabatic 
contraction to the dark matter improves the fit slightly and con- 
sistently yields a lower M^ILr ratio. Based on the range of de- 
rived M, ILr ratios and comparison to a range of literature IMFs, 
we argue that low M^ILr fits are preferred on physical grounds, 
implying that adiabatic contraction has operated in this system. 
These conclusions are robust to the presence of a possible central 
X-ray point source. 

Clearly the most significant source of uncertainty in our anal- 
ysis is the IMF. Observation of more extreme fossil systems may 
allow unambiguous detection of the adiabatic contraction effect 
on the dark matter, and provide evidence that these systems are 
indeed older than normal systems. As it stands, in the presence of 
such excellent data, better observational and theoretical progress 
on the IMF is necessary to draw definitive conclusions. 
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Appendix A: Background processing for 
RXC J231 5.7-0222 

In the case of RXC J23 15.7-0222 the source emission fills the 
XMM-Newton field of view, limiting the application of our stan- 
dard method, which relies on a source-free region from which 
to estimate the local background, and ultimately, the CXB con- 
tribution to the annular cluster spectra. For this system, we thus 
adopted the following procedure to calculate the surface bright- 
ness and temperature profiles. 
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A.1. Surface brightness profile 

We take the surface brightness profile subtracted from the FWC 
background data a nd fit it with the analytical model proposed 
by Vikhlinin et al. ( |2006|) plus a constant to estimate the CXB 
background level. We then take the mean between the CXB esti- 
mated with our standard method (which is an overestimate) and 
the constant level estimated from the fit. We then checked our 
estimate of the CXB contribution using ROSAT data from a re- 
gion around the position of RXC J23 15.7-0222, finding that the 
estimates are consistent. We use the mean CXB estimate and add 
the difference between the result from the standard method and 
the analytical model plus constant model fit as a systematic error 

A.2. Spectrai analysis 

An estimate of the spectral contribution from the CXB was ob- 
tained by simultaneously fitting the two outermost annuli with a 
model consisting of two unabsorbed thermal components corre- 
sponding to the galaxy and the local bubble, an absorbed power- 
law with slope fixed to F = 1 .4 corresponding to the contribution 
from unresolved AGN, and an absorbed thermal emission model 
representing the group component. Physical background param- 
eters were linked between the two areas (temperatures, power- 
law index, etc) but the normalisation of the background contribu- 
tion was allowed to vary. The final best fitting background model 
is then used in the fit to the inner annuli, with its normalisation 
scaled to the ratio of the areas under consideration. 

We checked that the CXB contribution estimated from the 
spectral analysis was consistent with that obtained from the sur- 
face brightness profile. 

Appendix B: Estimating the AGN contribution for 
RXC J021 6.7-4749 

RXCJ0216.7-4749 was found to have an X-ray bright AGN at 
its centre, manifesting itself after our standard deconvolution in 
a very cusped density profile and a flat inner temperature pro- 
file (see Fig. IB. 11 1. Various tests confirmed our suspicion; for in- 
stance an F-test showed that the addition of a power law com- 
ponent in the central regions considerably improved the spectral 
fit, and the surface brightness pr ofile was better fitted with an 
AB model (P ratt & Arnaudl2002h plus a point source convolved 
with the PSF 

We first sampled the central region in such a way as to max- 
imise the AGN contribution to the innermost temperature profile 
bin. In this 9" region, a single thermal emission model fit yields 
a;t'^/dof~ 3 and a temperature of kT ~ 2.5keV. Adding an ab- 
sorbed powerlaw improves the fit considerably (confirmed with 
an F-test), and suggests that the AGN contributes 40-50 per cent 
of the total counts in this region. We used several different meth- 
ods to estimate the AGN contribution, as listed below. 

1 . Fitting the surface brightn ess profile with an analytical 
model dVikhhnin et alj|2006l) plus a point source: 

- fit the surface brightness profile with the model (fossil 
plus point source); 

- integrate each component of the result (fossil and AGN) 
in each annulus; 

- fit the spectra of the annuli and tune the powerlaw nor- 
malisation so as to find a consistent count rate for the 
powerlaw component in the 0.3-2 keV band (correspond- 
ing to the extraction band of the surface brightness pro- 
file). 




100 
Radius (kpc) 

Fig.B.2. The mass profile of RXC J02 16. 7 -4749. Data points 
are the mass profile derived from hydrostatic equilibrium and 
solid lines are the best fitting model in each case. Black: stan- 
dard analysis with no correction for the AGN contribution (the 
inner point is omitted because the single thermal model gives 
a x^/y > 2). Green: analysis after excision of the inner 18" 
of mass profile data. Red: analysis with AGN contribution cor- 
rected as described in this Appendix. 



2. Fitting the surface brightnes s profile with an AB profile 
model (iPratt & Arn aud 2002) plus point source, repeating 
the steps in 1 above. Thi s alternative was introd uced because 
the analytical model of I Vikhlinin et al.l (l2006h is designed 
for cool core systems, whose brightness distributions behave 
differently to a point source when convolved with a PSF. 

3. Fitting the spectrum of the central region with a thermal 
emission model plus an absorbed powerlaw: 

- fit the spectrum of the first annulus with a source plus 
powerlaw model. The free parameters are the fossil prop- 
erties (temperature, abundance,...) and the normalisation 
of the powerlaw component. The powerlaw index is fixed 
to r = 1 .4. We use as a contraint the fit result of the first 
annulus only because the AGN contribution decreases 
rapidly outside this region and the spectral fit tends to 
find zero contribution at odds with the expectation from 
the surface brightness analysis; 

- estimate the AGN count rate in the 0.3-2 keV band from 
the spectral fit; 

- fit the surface brightne ss profile with an analytical model 
JPratt & Arnaudll2002l) plus point source, with the point 
source normalisation tuned to obtain an AGN contribu- 
tion consistent with the results from the spectral fit. 



We consider the final method to give the best estimate of the 
AGN count rate contribution as it combines both imaging and 
spectroscopic constraints. Using this method the AGN contribu- 
tion in the first three temperature profile annuli of RXC J02 16.7- 
4749 is estimated to be 40, 18 and 6 per cent, respectively. To 
compute the uncertainty on the AGN contribution we change its 
normalisation in XSPEC so as to obtain Ax' ~ 1 . fixing all the 
other parameters. The resulting relative error is ~ 3.5 per cent. 
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Fig.B.l. Surface brightness (left), density (middle) and temperature (right) profles of RXCJ0216. 7-4749. Black: standard analysis 
with no correction for the AGN contribution (the inner point of the temperature profile is omitted because the single thermal model 
gives ax^/v > 2). Green analysis after excision of the inner 18" of profile data. Red: analysis with AGN contribution corrected as 
described in this Appendix. 



Appendix C: Effect of central regions on mass 
profile modelling of RXCJ021 6.7-4749 

In this Appendix, we are interested in the effect of various 
ways of dealing with the AGN on the parameters resulting from 
an NFW model fit to the mass profile. The three cases we 
considered are illustrated in Figure IB. II which shows the sur- 
face brightness (left panel), deconvolved density (middle panel) 
and deconvolved temperature (right panel) profile of each case. 
Figure IB. 21 shows the resulting mass profiles and best fitting 
NFW models; model parameters are listed in Table IC.ll The 
three cases are: 

- No correction (black profiles). Here the point source is 
clearly visible in the surface brightness profile, and our stan- 
dard deconvolution process yields a strongly peaked density 
profile. The temperature profile is flat in the inner regions. 
The resulting mass profile has an excess of mass in the cen- 
tre compared to the NFW model fit, which does not fit the 
data points very well Ot'^/dof = 15.0/12). This fit gives the 
lowest value of M500 and the highest concentration. 

- Exclusion of the central 18" (green profiles). In this case we 
lose all constraints on the core mass profile of the group, 
relying on the outer data points to constrain the shape of the 
mass profile. The NFW model fit is good (y^/dof - 2.4/6), 
and yields values of M500 and c that are intermediate between 
a fit of the uncorrected profile and a fit of the profile with 
correction for AGN contamination. 

- Modelling of the AGN contribution as described above in 
Appendix |B] (red profiles). Here we have the tightest con- 
straints on the mass profile shape, from a profile that extends 
from [0.01 - 1] R500. The NFW fit is again good (y^/dof = 
10.0/13), but yields the highest value of M500 and the lowest 
value of c. 



NFW fit parameter TOT CUTCENTRE AGNMOD 

M500 (lO'^MO) To™ "~^"^ -^^- 



C500 



T c 7+0.39 



1.16+"" 



1.3 



;+().li 



2.22 



-(l,U 

-.+(17' 



+'0^5 



Table C.l. Eff^ect of excising the central regions (CUTCENTRE) 
or modelling the AGN contribution to the central regions 
(AGNMOD) on the parameters derived from an NFW model fit 
to the mass profile of RXC J0216.7-4749. 



ture profiles were thus corrected for the presence of this possi- 
ble central point source as described in Appendix IBJ The result- 
ing mass profiles were refitted with the mass models described 
in Section 15.21 The best fitting model parameters are given in 
TableEI] 
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Appendix D: Mass modelling results for 

RXC J231 5.7-0222 with correction for additional 
central point source 



As discussed in Section 16.1.41 a fit to the spectrum of the cen- 
tral region suggests the possible presence of a point source con- 
tribution at the < 10 per cent level. The density and tempera- 



16 



J. Democles et al.: Adiabadc contraction in fossils 



Table D.l. Results of fits to the mass profile of RXC J23 15.7-0222 with correction for a possible central point point source. 



Model 



C500 



M500 (10"Mo) i?5oo (kpc) 



/ga.;.500 



MJLr 



M5Qo/M',,'„ 



X /dof dvin, 



V 


7.14 ±0.60 




3.80 ±0.25 


V+star 


6.63 ± 0.61 




4.18 ±0.13 


V+ k star 


5.71 ±0.81 




4.18 ±0.13 


VAC+star 


5.31 ±0.29 




4.18 ±0.13 


V*AC+ k star 


5.70 ±0.14 




4.18 ±0.12 




5.27 ± 0.84 


0.118 ±0.021 


4.40 ±0.31 


+star 


5.43 ± 0.92 


0.161 ±0.037 


4.39 ± 0.23 


+<:star 


5.60 ±1.07 


0.177 ±0.092 


4.35 ± 0.22 


*AC+star 


4.88 ± 0.66 


0.262 ± 0.063 


4.25 ±0.21 


*AC+^ star 


5.42 ± 1.14 


0.188 ±0.104 


4.38 ±0.19 



509 ± 10 0.089 : 

523 ± 5 0.086 : 

524 ± 5 0.086 : 

523 ± 5 0.086 : 

524 ± 5 0.086 : 



0.011 
0.008 
0.008 
0.008 
0.008 



532 ±13 0.084 ±0.011 

532 ± 10 0.084 ± 0.010 

530 ± 9 0.084 ± 0.009 

526 ± 9 0.085 ± 0.009 

532 ± 8 0.084 ± 0.009 



0.72 ±0.07 18.3/13 0.55 

1.84 0.78 ±0.06 9.7/13 0.28 

3.34 ±1.04 0.78 ±0.06 7.7/12 0.25 

1.84 0.78 ±0.06 9.29/13 0.36 

1.38 ±0.31 0.78 ±0.06 7.34/12 0.25 

0.82 ±0.08 6.48/12 0.18 

1.84 0.82 ±0.07 6.57/12 0.20 

1.58±1.97 0.81 ±0.07 6.51/11 0.21 

1.84 0.79 ±0.07 7.31/12 0.30 

0.80 ±0.94 0.82 ±0.07 6.36/11 0.21 



